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Reflector Oriented Amplitude Recovery (ROAR) 

Vetle Vmje, NORSAR, 
Sist modifisert: May 3, 2000 4:14 pm 

1 Innledning 

1 seismikk er malet a kartlegge fysiske egenskaper i undergrunnen. Seismisk energi i form av 
lydb0lger blir sendt ned i undergrunnen og responsen (refleksjoner, diffraksjoner, transmisjoner) 
blir registrert av sensorer som oftest plasseres p5 overflaten av jorda, men som ogsS kan plassercs 
pa havbuiinen eller i borehull. 

Etter at signalene fra sensorene har blitt "prosesserf' (bearbeidet i datamaskin), kan geofysikere 
og geologer med dette som utgangspunkt lage karter over strukturer i imdergrunnen. Disse 
stnikturene representeres fcnrst og fremst av grenseflater mellom bergarter med forskjellige fysiske 
egenskaper, og disse grenseflatene blir registrert og kartlagt fordi de reflekterer seismiske balger. 

Denne patententseknaden beskriver en ny metode der en ved hjelp av seismiske mSlinger kan 
beregne refleksjonsstyrken, uttrykt ved den vinkelavhengige refleksjonskoeffisienten, ved disse 
grenseflatene. For a imnga misforstaelser vil de fleste faguttrykk skrives pa engelsk, i kursiv. En 
full forstSelse av denne patents0knade!n fonitsetter en viss kjennskap til geofysikk generelt og 
seismikk spesielt. 

2 Kjent teknikk pa omradet 

Det finnes en omfattende litteratur pa metoder som gar ut pa a finne refleksjonsstyrken til 
reflektorer i undergrunnen. En fellesbetegnelse som ofte blir brukt om disse metodene er "True- 
Amplitude Migration" (TAM). Det er ikke mulig innen et dokument som dette a gi et utfyllende 
resyme av alle metoder og varianter av metoder som er publisert, men et grunnleggende trekk ved 
de eksisterende metodene er at de alle S0ker a finne bade grenseflatens jposisjon og dens 
refleksjonsstyrke i en og samme operasjon. 

En oversikt over og sammenligning av noen eksisterende metoder kan finnes i [1]. Sentrale 
publikasjoner innen TAM er [2], [3], [4], [5], [6] og [7]. 

For a beskrive den nye metoden som det S0kes patent for her, ma vi begynne med en beskrivelse 
av Kirchhoff prestack dypmigrasjon ([5], [6], [7]), som ofte forkortes til PSDM (Pre-Stack Depth 
Migration). Vi tar utgangspunkt i Figur 1 som viser en ukjent reflektor, et generelt varierende 
hastighetsfelt v(x,y,z) og et skudd-mottakerpar med tilh0rende observert tids-trase. Et seismisk 
survey inneholder ofte millioner av slike skudd-mottakerpar. Merk at vi betrakter et 3- 
dimensjonalt medium (undergrunn), selv om illustrasjonene for enkelthets skyld er 2- 
dimensjonale. Vi antar na at PSDM blir brukt for a flytte (migrere) registrerte data fra 
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registreringene pS overflaten ned til en posisjon, x , i undergrunnen. En slik dypmigrasjon er 
generelt gitt ved formelen 

der / er antall traser som adderes i en generell posisjon x i undeigrunnen, 
tij{x) er gangtid #/ av de mulige stralebaner fta Skudd via x til Mottaker for 

skudd/mpttakerpar # z, 
p-[t-jix)] er den tidsderiverte av den /'te registrerte trase for gangtid ^y(x) , og 

w^j[x] er vekten som legges pa traseverdien Pflt^ix)] . 

De ukjente sterrelsene pa heyre side av Ligning (1), f,y(x) , J(i) og w,y[x] , blir beregnet ved 
hjelp av stralemodellering ([8], [9], [10]) eller endelig-differeansemetoder ([13], [14]). 

Nar et stort antall registrerte traser blir addert som beskrevet i Ligning (1), blir de reflekterte 
pulsene fra den ukjente reflektoren plassert ut i dypet langs isokronflater som alle tangerer den 
ukjente reflektoren. Dette vil f0re til konstruktiv intereferens langs reflektoren, og destrutiv 
interferens vekk fra reflektoren. Derved vil den migrerte seismikken, Fq^^{x), \isQ h0yere 
amplityder pa, og i naerheten av reflektoren slik at dens posisjon kan bestemmes. 

1 mye av utviklingen rundt PSDM og TAM ([2], [5], [6], [7]). er malet a finne vekter Wy[x] som 
gir amplitydeverdier /'obs(^) proporsjonale med refleksjonskoeffisienten (RK.) langs 

reflektorene. Refleksjonskoeffisienten [1 1] er forholdet mellom reflektert amplityde og 
innkommende amplityde til en plan b0lge som reflekteres fta en plan grenseflate. 

1 tillegg bar det vaert lagt mye arbeid i a finne ut hvilken av de J(i) mulige stralebaner fra Skudd 
via x til Mottaker som gir best resultat i migrasjonen. Et alment akseptert og brukt kriterium for 
et slikt valg er S brake den stralen fra skudd/mottakerpunkt til dyp-punkt x som bar sterst 
amplityde [1 2]. 1 det siste bar ogsa stralen med kortest gangvei vasrt valgt med gode resultater. 
Ved a velge kun en av de 7(0 mulige stralebanene fra Skudd via x til mottaker kan Ligning (1) 
omskrives til 

/ 

der t-{x) er gangtiden til den utvalgte stralebanen fra Skudd via x til Mottaker for 

skudd/mottakerpar # i , 
PiUjW] er den tidsderiverte av den fte registrerte trase for gangtid ^,(x) og 
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w,[x] er vekten som legges pa traseverdien /?/[f,(x)] . 



1 beskrivelsen av den nye metoden nedenfor, antar vi at vekten iv,[x] er 1 .0 for alle x og i. 

3 Beskrivelse av den nye metoden 
3-1 Oversikt over fbrkortelser 

For & lette lesingen av patentbeskrivelsen, folger en liste over forkortelser: 



3.2 Sammendrag og essens 

Metoden som det sekes patent for kalles "Reflector Oriented Amplitude Recovery" (ROAR). 
Ut g an gspunkt for ROAR er 

• Et seismisk survey^ dvs posisjonene til alle skudd og mottakere i en seismisk unders0kelse. 

• Karakteristiske egenskaper til skudd og mottakere (kildepuls, retningsavhengighet, 
filteregenskaper etc.) 

• Fre-stack dypmigrerte constant-offset eller constant-angle seismiske datakuber. 

• Posisjon og form til '^Target Reflector^\ dvs flaten der refleksjonskoeffisienten skal finnes, 

• Modcllparametre (seismiske hastigheter, tettheter, dempningsfaktorer, grenseflater etc.) til 
den delen av modellen som ligger over horisonten. 

Resultatet ved bruk av ROAR er 

• Diskret estimat av den vinkel-avhengige refleksjonskoeffisienten (RK'en) for alle belyste 
omrader av flaten (Ligning (25) og Ligning (31)). 

• En vektfunksjon (Ligning (47)) definert rundt alle diskrete verdier av RK'en med tilh0rende 
vinkler (Ligning (48)). 

Fordelene ved bruk av ROAR i forhold til eksisterende metoder er: 

• En handtering av reflektorenes krumning slik at bedre estimater av RK'er oppnas. 

• Diskrete estimater av RK for belyste omrader pa flaten med tilh0rende refleksjonsvinkler. 
Dette gir samtidig en kvantitativ oversikt over darlig og manglende belysning av flatene ved 
alle refleksjonsvinkler. 

' • En vektfunksjon for alle RK estimatene som gir oppl0sningen langs Target Reflector til hvert 
enkelt estimat. 

• Mulighet for a finne RK'ene i modeller med bade akustiske og elastiske modcllparametre. 

• Mulighet til a inkludere forskjellige b0lgemoder for a estimere RK'ene for PP-, PS-, SP- og 
SS-refleksjoner. 

• Muligheter for a korrigere for forskjellige typer av multipler ved a inkludere disse i 
stralemodelleringen. 



RK 
ROAR 
PSDM 
TAM 



Refleksj onsKoeSisient 

Reflector Orienter Amplitude Recovery 

Pre Stack Depth Migration 

True Amplitude Migration 
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• En bedre handtering av begrenset utstrekning av det seismiske^wrvey og huUer og mangier i 
surveyet. 

• En generell mulighet til i forbedre estimatet av RK'en ved a inkludere mer komplett 
stralemodellering som f.eks. kan baseres anisotrop straleteori og anelastisk dempning. 

• Redusert beregningstid og mindre bruk av diskplass. 

Som nevnt i kapittel 2 seker de eksisterende metodene 4 l0se dette problemet ved S finne bade 
flatenes posisjon og dens RK i en enkelt prosess, nemlig ved § bruke forskjellige varianter av 
PSDM med vekter pa traseverdiene. I motestning til dette finner ROAR posisjon og RK ved en 
prosess i fern trinn. Trinn 1, 2 og 3 er kjent teknikk, mens Trinn 4 og 5 er essensen i oppfinnelsen: 

Trinn 1: De registrerte trasene migreres vhja PSDM til dyp etter Ligning (1). Det er en fordel, 

men ikke nodvendig, 5 bruke den enkleste og raskeste versjonen av PSDM, der vekten 
er konstant 1 .0. Resultatet av denne forste prosessen er en dyp-migrert seismisk kube, 

^obs(^)- 

Trinn 2: /*obs(^) ^o^^^^ civs at posisjonene til reflektorene i undergrunnen blir funnet. 

Reflektorene representeres som flater som markerer diskontinuerlige sprang i de 
fysiske egenskapene til bergartene. Utifra disse flatene og de seismiske hastighetene 
(de samme hastlgheter som ble brukt i PSDM prosessen i Trinn 1) bygges en 
jordmodell opp i datamaskinen. En av flatene i modellen velges ut som en "malflate" 
( Target Reflectof), 

Trinn 3: Ved hjelp av stralebasert modellering beregnes de reflekterte pulsene fra Target 

Reflector for alle skudd/mottakerpar i hele surveyet. Siden vi simulerer seismikken i en 
jordmodell i datamaskinen, kan vi fritt sette RK'en til Target Reflector \\\ en konstant 
verdi, f.eks 1 .0. Nar RK'en settes til 1 .0 vil all energi i pulsen reflekteres fra Target 
Reflector 

Trinn 4: Vi kjorer en lokal PSDM pa de syntetiskd trasene fra den stralebaserte modelleringen i 
Trinn 3 i et band rundt Target Reflector sWk at vi far et modellert PSDM resultat 
^Mod(^) ^^^^^ Target Reflector. Vekten (se Ligning (1)) som brukes i denne lokale 
PSDM ma i hovedsak vaere lik vekten i Trinn 1 . 

Trinn 5: Vi foretar malinger av amplitydene ved f.eks. "kvadrat" eller "norm" (se kapittel 3.5 og 
3.6) fra den reelie, dypmigrerte PSDM kube /'obs(^) dividerer med de tilsvarende 
malingene fra modellert PSDM PuodW ■ Det kan vises at resultatet er et estimat av 
den vinkelavhengige RK. Dette gjentas for punkter langs hele Target Reflector, og et 
kart over RK for Target Reflector kan lages. 

I den detaljerte beskrivelsen nedenunder, settes det i kapittel 3.3 opp et matematisk uttrykk for 
PSDM av observert seismikk, dvs for den dypmigrerte seismiske kuben,Pobs(^) • tilsvarende 
uttrykk for PSDM av den modellerte seismikken, Pp4o(j(x) . utledes i kapittel 3.4. 1 de pafalgende 
kapitler beskrives to varianter for a finne den vinkelavhengige RK vedTarget Reflector samt 
formlene for a finne tilherende vektfunksjonen (oppLasning) og refleksjonsvinkel. 
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3.3 Uttrykk for observert seismikk 

Vi skal na finne et altemativt uttrykk for Pq^^O^) , resultatet av en PSDM fra Ligning (2) i 
omradet rundt en utvalgt reflektor i modellen. Uttrykket gir /'obs(^) ™ 
straleteoretiske komponenter. 

Anta at ^obst^) ftonet for alle x , ved Ligning (2) for et bestemt q^e^-intervall, dvs for et gitt 
utvalg av skudd/mottakerpar som inng^r i sunimasjonen i Ligning (2). Pq^^O^) er da en constant- 
offset seismisk kube. Den seismiske kuben blir tolket, og eksplisitte reflektorer i form av romlige 
flater blir lagt inn i modellen. En av disse reflektorene velges ut som en Target Reflector, 

Definisjon: En reflektor der en 0nsker k estimere refleksjonskoeffesienten (RK) vhja ROAR kalles en 
Target Reflector 

Det antas at stralemodellering ([8], [9], [10]) i en jordmodell som inkluderer Target Reflector kan 
brukes til k estimere trasene, Pi[t]/\ Ligning (2). Ved stralemodellering kan man beregne pulsen i 
mottakeren fra en bestemt strale # k som er reflektert fra Target Reflector med falgende ligning: 

g^[t] = K^{RkQk)s{t-x^)-ha{R^Q^}hit-Xf^) (3) 
der Xj^. er gangtiden til stralen, 

s{t) and h{t) er kildepulsen og dens Hilbert transformasjon, 
5(0 1 ^0 [0, r] 



1 ^0 ti 

► J =0 n 



, T et kildepulsens varighet i tid, 

h{t)\ =0 /€£ [0, T] 

R^ = /?^((pjt) erdenukjente,kompleksevinkeIavhengigeRKtilstralenirefleksjonspunktetR 
pa Target Reflector 

og Qf^ er et komplekst tall som gir den samlede effekten av geometrisk spredning, 

transmisjons- eller refleksjonstap ved andre grenseflater (etm Target Reflector)^ 
faseskift ved at strSlen gar gjennom fokalpunkter/linjer {caustics) og anelastisk 
dempning. 

Ved stralemodellering kan en rekke parametre beregnes langs strSlebanen, 1 tillegg til 
refleksjonspunktet, R, kan stralens irmkomstretning til R, utgangsretning fra R, utgangsretning fra 
Skudd, innkomstretning til Mottaker og en lang rekke andre parametre beregnes, Bruken av disse 
parametrene blir diskutert nedenfor. 

Som nevnt ovenfor, kan hver av de / enkelt-trasene som adderes inn i x i Ligning (2) estimeres 
ved stralemodellering 

PiUiW] = J^g.lhW] (4) 

k= 1 
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der P(i) er antall strSler som er fuiinet mellom skudd/mottakerpar nummer z. 

Vi setter sa Ligning (4) inn i Ligning (2), og antar at vekten w -[x] er 1 .0, og vi fir 

^Obs(X) = S E S,UiW] (5) 

Vi vet at det har vasrt en konstruktiv interferens langs Target Reflector mellom de migrerte, 
registrerte pulsene. Det betyr at det ikke er nedvendig a addere traser for verdier av x som ligger 
utenfor en viss avstand fra Target Reflector, Dette vil rednsere antall pulser, g^, som brukes i 
. Ligning (5). Dersom vi antar at X ligger innenfor et band rundt Target Reflector som vist i Figur 
2, kan Ligning (5) uttrykkes som 

N 

^Obs(x) = ^gniW^ xe X^ (6) 

der N er antall strale-beregnede enkeltpulser innenfor bandet X g x j- som er slik at 

'C„<^„(x)<T„ + r xexy. (7) 

der T er varigheten til pulsen, som vist i Fignr 2. Fra denne figuren kan en se at kravet i Ligning 
(7) er oppfylt for plasseringen av x i Figur 2 fordi x i dette tilfelle ligger mellom de to 
isokronene. 

Vi introduserer sa en ny definisjon, en Stacking Stick, 

Definisjon: En Stacking Stick (# J) er en kurve med en bestemt lengde som krysser Target Reflector 
som illustrert i Figur 3. Den har en bestemt lengde. som gir bredden p5 bSndet rundt 
Target Reflector. veiges stor nok til S dekke noen bolgelengder av Pq^^O^) rundt 
Target Reflector. En Stacking Stick er parametrisert vhja Xy(^y) som markerer den 
romlige posisjonen. Argumentet starter med 4y = ^ja i den ene enden av Sfac/f/ng 
Stick, passerer = 0 i krysningspunktet mellom Stacking Stick og Target Reflector, og 
slutter i = i den andre enden av Stacking Stick, 

Anta at x i Ligning (6) ligger pa en bestemt Stacking Stick ^ji^j) - ^ ^'^EM^ 3 er Stacking Stick #j 
tegnet som et rett, vertikalt linjestykke som krysser Target Reflector, men kan ogsa f eks plasseres 
perpendikulasrt til Ta/^er^^ecror. /'obs(^) eroftest representert pa et regulaert, 3-dimensjonalt 
nettverk (i.e. en seismisk kube) der de diskrete x-posisjonene er gitt ved 

^ijk = X^+<(z-l)c/jc,0-l)dy,(A:-l)rfz> (8) 
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der er posisjonen til det ene hjemet av nettverket og {dx, dy,dz) er avstandene mellom de 
diskrete x-verdiene i hhv x-, y- og z-retning. For a finne verdien av Pq^s(x) i stacking stick #7, 
brakes interpolasjon (lineser elle h0yere ordens) slik at vi kan fiiine Pq^^{Xj(^j)) . 

For enkelthets skyld betegner vi den observerte trasen i stacking stick # j med 

^Obs(^y) = ^Obs(X/^y)) W 

Samme forenkling kan brakes i betegnelsen av Kirchhoff gangtiden Skudd->x ->Mottaker, nar 
X = Xy(^y) liggerpa Stacking Stick if J: 

Vi introduserer sa f^i^j) , en tilnaermelse til Kirchhoff gangtid i Ligning (6). Tilnsermelsen 

t\i^j) bruker sakalt paraksial ekstrapolasjon [8] av gangtidsverdiene rundt R (se Figur 3) for 
stralebane # « fra Skudd, via x^.(^y) pa Stacking Stick j til Mottaker. Paraksial ekstrapolasjon 
innebaererat 

t\aj)-tni^j) (11) 

for Stacking Sticks i nserheten av R. Ved bruk av denne tilnaermelsen kan fra Ligning 

(3) og Ligniag (6) uttrykkes som 

gj^(^y)] = R^{RMQ,}sit\(^j) - x„) - lm{R^(i?^)QJh(t\aj) - T„) (12) 

Siden refleksjonsvinkelen, <p„ , inngar i Ligning (12), kan vi introdusere vinkelavhengigheten i 
^Obs(^y) omskrive Ligning (6) til 

Pobs(^y.9y) = X [Re{i?„(9„)e„}i(^{^y)-X„)-Im{/2„((p„)0„}A(/'„(4y)-X„)] (13) 
rt = 1 

der (pj er et veiet snitt av de N refleksjonsvinklene (p„ i Ligning (13), Vi kommer inn pa en 
naermere definisjonen av denne vinkelen senere. 

Vi antar en reell RK, i?„(q>„) , som f.eks ved PP refleksjon ved underkritisk vinkel. Vi utnytter 
denne egenskapen til R^((p^) i Ligning (13) og far 

N 

^Obs(^,>^y)= J,KWPn(^j) (14) 

der 

Pni^j) - Re{0Ji(^',(^,-)-x„)-lm{e„}A(r'„(4y)-T„) (15) 
Vi introduserer sa en gjennomsnittlig RK, i?(^y, (py) : 



Side 7 av 18 



PATENTS0KNAD: Reflector Oriented Amplitude Recovery (ROAR) 



N 



n= 1 

der R(^j, cpy) kan fmnes ved S kombinere Ligning (14) med Ligning (16): * 

N 

^(^y ^y) = X ^«(^y)^n(9„) (17) 
«= 1 

medvekter, vv„(^y), 

1 = 1 

3.4 Uttrykk for modeliert seismikk 

Vi har na fiinnet en stralebasert tilnaermelse til PSDM av den observerte (eller registrerte) 
seismikken gitt i Ligning (16) for et band med bredde langs Target Reflector, Wi skal na finne 
uttrykket for (strale) modellerte traser i Stacking Stick #j i en modell der Target Reflector har en 
RK som er definert til 1 .0 for alle vinkler o for alle punkter. men som ellers er lik den virkelige 
modellen. Ved samme resonnement som ved utledningen av Ligning (6) Sr vi: 

n = I 

der t\(l,j) er Kirchhoff gangtiden beregnet ved paraksial ekstrapolasjon [8] rundt R (se Figur 
3) for stralebane # « fra Skudd, via y^ji^j) pa Stacking Stick # j til Mottaker , 

g'„[f„i^j)] er den tidsderiverte av den «'te stralemodellerte trase g'„[r'„(^y)] som uttrykkes 
ved a sette R^ = 1.0 i Ligning (3): 

Ved a sammenligne Ligning (20) med Ligning (15) ser vi at 

som settes inn i Ligning (19) slik at vi far 
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som er uttrykket for det modellerte resultatet av u-vektet PSDM i en modell der RK'en settes lik 
1 .0 p5 Target Reflector. Antall straler N er samme mate som i Ligning (7) ovenfor, utvalget av 
alle straler mellom alle skudd/mottakerpar som er slik at 

for minst ett av punktene ^ji^j) Stacking Stick #7. 

3.5 Estimering av refleksjonskoeffisient, Kvadratmetode 

Utfra forutsetningene gitt for utledningen av Ligning (16) og Ligning (22), kan det vises at et 
estimat av kvadratet av den vinkelavhengige RK'en ved krysningspunktet mellom Stacking Stick 
# j og Target Reflector er gitt ved forholdet 



= 0, {pj)f = ; (24) 



der 9l£:(^y = 0, (py) er estimatet av RK'en for krysningspunktet mellom Stacking Stick #j og 
Target Reflector for refleksjonvinkel <py, ^obsC^y ^y) observerte trasen i Stacking Stick #j\ 

^Mod(^y) er den stralemodellerte trasen i Stacking Stick it j\ F(^y) er en valgt vektfiinksjon langs 
Stacking Stick #j og (^y^, ^j^) er ekstremalverdiene til langs Stacking Stick #J (se Figur 3). 

Den diskrete versjonen av Ligning (24) er som f0lger: 



-Mod 

der 

L 



[91^(^-0,9/)]' = 1^ (25) 



^Mod = E {^Mod(^y;t)^(^y)t)}'A^ (26) 
k= 1 

Og 

L 

^Obs = Z {^obs(^y)t.^y)^(^yit)}'A^ (27) 
k= 1 

der er en diskret verdi av , ^^ji^jk) tilhenrende posisjon # A: til Stacking Stick #j og 

er avstanden mellom l^ji^jk), \'(^jik+ i ))' - ^^^S>^ Stacking Stick #y. 
Med Ligning (16) og Ligning (22) innsatt i hhv Ligning (26) og Ligning (27) kan det ogsa utledes 
en vektfunksjon for bidraget fra alle de individuelle RK'ene rundt Stacking Stick #j\ 
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Det kvadratiske estimatet til RK'en kan uttrykkes som 



N 



(28) 



der vektenpS RK'en, /?„(9„X i refleksjonspxmktet, R, til strale # n er gitt ved 




(29) 



^Mod ^ Ligning (26) og den dypmigrerte puis PjXjk) er gitt ved resultatene av 

stralemodellering og Ligning (15). P^i^jk) beregnes fra en strale med refleksjonsvinkel (p„ som 
vistiFigurS, 

3.6 Estimering av refleksjonskoeffisient, Norm-metode 

Utfra forutsetningene gitt for utledningen av Ligning (1 6) og Ligning (22), kan det vises at et 
estimat av absoluttverdien av den vinkelavhengige RK'en ved krysningspunktet mellom Stacking 
Stick # j og Target Reflector er gitt ved forholdet 



der 9?jv(^y = 0, ^y) er estimatet av RK'en for krysningspunktet mellom Stacking Stick #j og 
TflA^ef/^e^ecror forrefleksjonvinkel (Py , ^obs^^y ^y) erdenobservertetraseniSracfa/i^S'ftc 
/, Z'ModC^y) den stralemodellerte trasen i Stacking Stick #7, F(^y) er en valgt vektfunksjoi 
langs Stacking Stick # J og {^j^, ^y^) er ekstremalverdiene til langs Stacking Stick #J (se 
Figur 3). 

Den diskrete versjonen av Ligning (30) er som f0lger: 




jJ'>Obs(^y.9y)/=-(4y)|d^y 



(30) 




Obs 



A'; 



(31) 



Mod 



der 



L 



N 



Mod 



= E I^Mod(M^(M^^ 



(32) 



Og 
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L 



^Obs = X |^Obs(^J^^9y)^(^y^)lA^ (33) 



k= I 



der er en diskret verdi av , ^^ji^jk) tilhcnrende posisjon # * til Stacking Stick #j og 
er avstanden mellom Xj{%jj^), x^C^y^j^^ i)), ... langs Stacking Stick #j. 

Med Ligning (16) og Ligning (22) innsatt i hhv Ligning (32) og Ligning (33) kan det ogsa utledes 
en vektfunksjon for bidraget fra alle de individuelle RK'ene rundt Stacking Stick #j. 

Absoluttverdien til estimatet av RK'en kan uttrykkes som 



N 



|9l^(^y = 0, 9y)| = i: IVni^j, ipj)RM (34) 
n = 1 

der vekten pa RK'en, R^((p„), i refieksjonspunktet, R, til strale # w er gitt ved 



^Mod Ligning (32) og den dypmigrerte puis Pni^jk) er gitt ved resultatene av 

stralemodellering og Ligning (15). P^i^jfc) beregnes fra en strale med refleksjonsvinkel (p„ som 
vist i Figur 3. 

3.7 Vekter og refleksjonsvinkler 

Vi har beskrevet to metoder innen ROAR for a estimere absoluttverdien til den vinkelavhenginge 
RK, , for et Target Reflector. De to er "Kvadratmetoden" (Ligning (25)) 

[9l£(4y = 0,9y)]' = 1^ (36) 

^Mod 



Og "Norm-metoden" (Ligning (31)) 



l^^i^j = 0, 9y)| = ^ (37) 
Mod 



- 2 

Det er vist ovenfor (Ligning (28) og Ligning (34)) at de estimerte RK'ene = 0, cpy)] og 

|SKy^(^y = 0, (py)| kau uttrykkes som veide gjennomsnitt av RK'ene i A^refieksjonspunkter i 
omradet nindt Stacking Stick # J. 

Det veide gjennomsnitt for a estimere RK'en i Stacking Stick # j for Kvadratmetoden er gitt ved 
Ligning (28): 
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n = l 

og tilsvarende for Norm-metoden er gitt ved Ligning (34): 

, N 

\^Ni^J = 0. (p,-)| - S f^^«(^y, (39) 
rt = I 

Formlene for vektene for hver av disse metodene er gitt i Ligning (29) og Ligning (35). Hvert av 
de N bidragene i Ligning (38) og Ligning (39) korresponderer til en bestemt strale med et 
refleksjonspunkt, R, og en RK i2„((p„) i naerheten av Stacking Stick # j\ AUe RK'ene i et omrade 
rundt Stacking Stick # j bidrar til estimatet 91 av RK'en slik at Ligning (38) og Ligning (39) er et 
godt utgangspunkt for a finne opplesningen (resolution) i RK'en. 

Vi antar na at Target Reflector er parametrisert ved et triangelnettverk som vist i Figur 4, 1 hver 
node i nettverket er det en Stacking Stick. Det antas at den virkelige RK'en er konstant i hver av de 
P trianglene som utgjcrr nettverket. 

T T T T 

Den (konstante) verdien til RK i hver av de P trianglene ct R^,R2,R^, .^>,Rp^ 

T T T T 

Til hver av RK'ene er det knyttet en gjennomsnittlig refleksjonsvinkel (pi , 92 ' ^3 ' 9/> 
definert i Ligning (42) nedenfor. 

Demest uttrykker vi Ligning (38) og Ligning (39) som 

/I = 1 



der 



= [9l£(^y = 0, (py)]^ og a„ - iVniJ^p (py) for Kvadratmetodcn 



= |9l^(^y = 0, (py)| og a„ = ^«(^y, (py) for Norm-metoden 

Hver av RK'ene /2„((p„) i ligningen korresponderer til et refleksjonspunkt i Figur 4. 
Sa definerer vi to undermengder: 
•' Sp : Undermengden av alle N refleksjonspunktene i Ligning (40) som ligger i triangel p , 

T 

• S : Undermengden av alle P triangler i Target Horizon som inneholder minst ett av de N 
refleksjonspunktene i Ligning (40) 

T 

Dermed kan vi finne relasjonen mellom i2„(q>„) og R^ : 

^„(q>J = < (41) 
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for alle ^„(9„) nied refleksjonspunkter i , i.e. for alle ns 

T 

assosiert til , 



T T 

Vi finner ogsa et estimat av refleksjonsvinkelen (p assosiert til RK'en 



cp^ = (42) 

#-elementeri 



Ved a bmke Ligning (41) i Ligning (40), kan Ligning (40) skrives som 



der 



= X (44) 

Vi 0nsker a finne en vektfunksjon som er "pr arealenhet" av Target Reflector. En kontinuerlig 
vektfiinksjon A(x^) ville matte integreres med en kontinuerlig RK fimksjon R (x^) pa fialgende 
mate 

b = jA{xs)R^{xs)dS Xs^ Target Reflector (45) 
5 

som kan diskretiseres ved bruk av var triangelparametrisering 

b = ^A^^AT^ (46) 

der ATp er arealet til triangel p og {A^, R^) er hhv vektfimksjonen og RK'en assosiert til 
triangel p. 

Sammenligner vi Ligning (46) med Ligning (43) ser vi at den diskrete tilnaermelse til en 
vektfiinksj on vil bli 



T 

for vekt per arealenhet for RJC'en i triangel p i undermengde S , 

* T 

Refleksjonsvinkelen (py assosiert til estimatet av b finnes ved et veid snitt av vinklene (p^ fra 
Ligning (42) 
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4 Figurer 




Isokron ^,-4(X) 
Isokron ^,-2(^) 



\^ Ukjent Reflektor 



Figur 1: 1 Kirchhoffprestackdypmigrasjon, adderes (stackes) den registrerte trase for skudd/mottaker par 
# / i dyp-posisjon X . / dette tilfelle er del 4 forskjellige strdlebaner som er kandidater til en reflektert puis 
registert i mottakeren: A-1 , A-2, B-1 . B-2. Vanligvis vil bare en av disse bli brukt i migrasjonen. Ifiguren er 
isokronene korresponderende til to av de mulige strdlebanene tegnei opp. Isokroner erflatersom svarer til 
konstant gangtid mellom skudd og mottaker. Der disse isokronene tangerer (den ukjente) reflektoren, har vi 
et refleksjonspunkt. I dette tilfelle er det to strdler tned gangtid hhv t^2W ^^(X) reflekterer i dyp~ 
punktet X. 
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Puis beregnet ved strSleteori 



< 



Skudd 



Isokron x» 




Isokron x„ + T 



Target Reflector 



Figur 2: 1 Kirchhoff prestack dypmigrasjon, adderes (stackes) den registrerte pulsen for skudd/mottaker 
par (U i) i alle punkter mellom de to isokronene, inkludert i dyp-posisjon X. Det antas at innenfor et band 
med en bestemt bredde L er det konstruktiv interferens mellom alle N strdlebaserte pulsene som legges ut. 
Derfor er det mulig a redusere antall pulser som migreres for et gitt punkt X innenfor bdndet. Bare de 
pulser med X innenfor de to isokronene vil stackes inn i X . 
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Skudd 




Registrert puis: g„it) 
Gangtid: ^ = x„ 



Refiektor 



Kirchhoff gangtid: - t„(^j) 

— Migrasjons-isokron, t\{^) 



Dypmigrert, modellert puis: 



Stacking Stick #7 



Target reflector 



Figur 3: Den modellerte puis # n for skudd/mottaker par # / adderes (stackes) i stacking stick # / En 
stacking stick er en kurve parametrisert med ^ - som hysser Target Reflector med parameterverdi 
t^ - — Q . i^'j^ er en diskret posisjon langs StacJdng Stick De andre stacking stick 'ene er plottet som korte, 
vertikale linjesegmenter langs Target Reflekton Stacking Stick 'ene avgrenser reflektorsegmenter med areal 



Konstant RK, konstant refleksjonsvinkel 




Refleksjonspunkter 



Figur 4: Target Reflector (sett i fiigleperspektiv) er parametrisert med et triangelnettverk RK'en i Stack- 
ing Stick # j er fldnnet utifra strdler med refleksjonspunkter som er markert med fylte sirkler. Den "virke- 
lige " RK'en pd flaten er parametrisert med konstant, vinkelayhengig RK, R^j hver av de P trianglene. 
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5 Utf0relsesformer/t^3ev^P^^'^ 

Se vedlegg "REFLECTOR ORIENTED AMPLITUDE RECOVERY". 

6 Referanseliste 

[1] Gray, S. H., 1997, True-amplitude seismic migration: A comparison of three approaches. 
Geophysics, 62, 929-936 

[2] Beyikin, G., 1985, Imaging of discontinuities in the inverse scattering problem by inver- 
sion of a causal generalized Radon Transform: J. Math. Phys,, 26, 99-108 

[3J Berkhout, A. J., 1 985, Seismic migration: Imaging of acoustic energy by wavefield 
extrapolation A. Theoretical aspects, Elsevier Science Publ 

[4] Berkhout, A. J., and Wapenaar, C. R A., 1993, A unified approach to acoustical reflec- 
tion imaging. Part 11: The inverse problem: J. Acoust Soc, Am., 93, 2017-2023 

[5J Hanitzsch, C, Schleicher, J„ and Hubral, P., 1 994, True-amplitude migration of 2-D 
synthetic data: Geophys. Prosp,. 42, 445-462 

16J Bleistein, N., 1987, On the imaging of reflectors in the earth. Geophysics, 52, 931-942 

[7] Hanitzsch, C, 1 995, Amplitude preserving prestack Kirchhoff depth migration/inversion 
in laterally inhomogeneous media: Ph,D. dissertation. University of Karlsruhe. 

[8] Cerveny, V,, 1 985, The application of ray tracing to the numerical modeling of seismic 
wavefields in complex structures: Handbookof Geophys. Expl, 15A, 1-124, Geophysical 
Press. 

[9J Vinje, V,, Iversen, E., and Gjoystdai, H., 1993, Traveltime and amplitude estimation 
using wavefront construction: Geophysics, VOL. 58, NO. 8; P.l 157-1 166 

[lOJ Vinje, V., Iversen, E., Aastebol, K., and Gjoystdai, H., 1996, Estimation of multivalued 
arrivals in 3D models using wavefront construction. Part 1 & 11: Geophysical Prospect- 
ing, 44, R8 19-858 

[11 J Aki, K. and Richards, P.G., 1980, Quantitative Seismoligy, Theory and Methods. Vol- 

imie 1, Chapter 5, ISBN 0- 71 67-1058-7 (v. 1 ) 
[12] Geoltrain, S. and Brae, J., 1993, Can we image complex structures with first-arrival 

traveltimes? Geophysics, 58, 564-575 
[13J Vidale, J. E., 1988, Finite-difference calculation of traveltimes: Bull Seis. Soc, Am., 78, 

no. 6, 2062-2076 

[14] Fodvin, P., and Lecomte, 1., 1991, Finite-difference computation of traveltimes in very 
contrasted velocity models: A massively parallel approach and its associated tools: Geo- 
phys. J. Int., 105,271-284 




Side 18 av 18 



REFLECTOR ORIENTED AMPLITUDE RECOVERY 



VETLE VINJE 

NORSAR, RO.Box 50, N-2027, Kjeller, NORWAY, vetle@norsar.no 



Summary 

I propose a new method for solving the classical inversion problem of finding the angle-dependent reflec- 
tion coefficients along selected target horizons in the subsurface. The input to the method is seismic data 
cubes from constant offset Kirchhoff migration and reflectors from interpretation on the cubes. One or more 
target reflectors are selected, and ray modelling by wavefront construction is done, creating ray data for all 
shot/receivers in the survey. Based on these modelling results, amplitude correction maps for different 
reflection angles are created and applied to the amplitudes measured on the seismic data. The corrected 
amplitudes are good approximations to the angle dependent reflection coefficients for all illuminated points 
on each target reflector. For each point, a weight function is computed which gives a quantitative resolution 
of the estimate of the reflection coefficient. 

Introduction 

In the exploration industry the ultimate goal is to obtain a precise and detailed knowledge of the lithology 
in the earth. An important source of such knowledge is the amplitudes from processed seismics which are 
assumed to be proportional to the (loosely termed) reflection strength of structures in the subsurface. This 
assumption is generally not true, and several authors have addressed the problem of finding the "true ampli- 
tudes" witiiin the last 20 years ([1], [2] and [3] to mention some). The metiiods described by these authors 
do not assimie any a priori knowledge of the reflectors since the amplitude corrections are incorporated 
direcfly into the migration methods. In the method described here, called "Reflector Oriented Amplitude 
Recovery", the migration and the estimation of the reflection coefficients (i.e., the inversion) are completely 
decoupled. This makes it possible to include the reflectors in the inversion and take their curvature and slope 
direcfly into account. 

Reflector Oriented Amplitude Recovery 

Consider Figure 1 which shows die recorded seismic trace for shot/receiver pair # /. A primary PP reflec- 
tion from a target horizon appears as a pulse on the seismic trace. In Kirchhoff depth migration, the trace 
(including the pulse) is distributed along migration isochrones which are surfaces with constant travel time 
shot->isochrone->receiver. The isochrone corresponding to the ray-theoretical traveltime of the pulse 
(shown in Figure 1) intersects one of the stacking sticks. A stacking stick is simply a portion of a depth trace 
wifli a length Lp and a parametrization x^(^) where 4 ^ [^^, f^] , tiiat crosses die target reflector. In 
Figure 1, the depth migrated recorded trace # i is added (stacked) into die stacking stick. This is repeated 
for all, or a subset (e.g. constant offset), recorded traces in die survey, giving die depth migrated trace cor- 
responding to a Kirchhoff migration. Eq. (1) gives the madiematical expression of the summed, unweighted 
pulse, P(|) , in the j'th stacking stick: 

M M 

^obs(^) = E ^/(^) = E 
i = 1 / = 1 

where r,-(^) is die traveltime of die Kirchhoff ray pafli (see Figure 1) in a position ^ in die stacking 
stick, PiUii^)] is die time derivative of die recorded trace p^ evaluated in travel time f .(^) and M is the 



number of traces contributing to the stack. 



1) 



Shot 





^- Traveltime 


Recorded trace: 






Receiver 




Reflector segment # n 

Stacking Stick # j 



Kirchhoff ray path 
Traveltinne: 

Migration isochrone 
migrated trace: = P,[^(^)] 

Target reflector 



Figure 1; The recorded trace for shot/receiver pair # i is stacked into stacking stick #j which is a part of a 
depth trace with a particular length. The other stacking sticks are plotted as ticks on the target reflector 
separating reflector segments of area AA^ . 

If we assume that the input pulses to this Kirchhoff migration are the result of ray tracing, the result, 
^ynth(^) ' *® stacking stick will be (usmg well-known ray theory, see e.g. [5]): 

i=\ i = 1 

where N is the total number of primary ray pulses stacked into the stacking stick, 
Tj. is the traveltime of the ray, 

s{t) and h{t) are the source pulse and its Hilbert transform, 

s{t) 1 ^ 0 r e [0, r] ^ duration of the source pulse), 

MO J =0 [0,r] 

R] is the (complex) angle dependent reflection coefficient of the ray at the reflection point, R, at the target 
horizon and Q/ is the (complex) combined effect of geometrical spreading and transmission loss of the ray 
from source to receiver (not to be confused by the geometrical spreading matrix). 

Within the approximation of ray theoiy, the true, stacked pulse, Pobs(^) . in Eq. (1) may be expressed in 
the same way as the synthetic pulse, ^^ynth(^) ' i^ obtain 
N N 

t = 1 / = 1 

where /?- = Ri{(^i) is the unknown angle dependent reflection coefficient in the real model in die reflec- 
tion point R (see Figure 1) of the z'th ray, and (p • is the reflection angle of the ray. ^ 

In order to find tfiis reflection coefficient we start by setting the syntfietic reflection coefficient in Eq. 
(2) to unity for all reflection angles. Then we compute the two quantities 



1 Jk= 1 

where are L discrete positions with spacing along the stacking stick. 

It can be shown that a good approximation to the angle-dependent reflection coefficient in the reflector at 
stacking stick j is given by 

SR^(cp') = 1^ (5) 

^ Synth ^ 

where the reflection angle (p' is a weighted sum of the N individual reflection angles (p • contributing to 
the estimate in Eq. (5). The ^yeight associated to the fth stacking stick for segment n, Wj^ , is computed 
using the input data cube, Pobs(^) ' *® individual ray-modelled, depth-migrated pulses P^(^) from 
Eq. (2). These weight functions give valuable information on the resolution of the estimate in Eq. (5), as 
will be shown in the example below. 



Synthetic example 

The simple 2D model used in the testing of the new concept is shown m Figure 2. 
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Figure 2: Two-layer model with variable P- and S-velociiy and density below target reflectonThe rays to 
every 4th receiver for shot # 700 of a total of 138 are shown. 



Migrated Depth Section. Offsets (25-300)m 
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Figure 3; Result of Kirchhoff migration of the 
lowest offset interval. 




The model consists of two layers separated by a tar- 
get reflector with a subduction. The lower layer has 
discontinuous changes in the seismic properties as 
indicated on the figure. 

At the surface of the model, a marine type survey of 
138 shots and 120 receivers per shot is modelled by 
wavefront construction [4]. The result is convolved 
with a 10 Hz zero phase Ricker wavelet to give the 
mput SEG-Y data to the test. Using these SEG-Y 
data, Kirchhoff migration is done in 10 offset inter- 
vals (offset bins) within die maximum offset of 3000 
m creating 10 seismic sections of which one is shown 
in Figure 3. Figure 4 shows the measured amplitude 
along die target reflector, which is obviously not a 
reliable way of finding die reflection coefficient. 
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Figure 4: Maximum amplitude in the seismic 
section in Figure 3 and true reflection coeffi- 
cient along the target reflector 



^^.-i-Estimated values, point p (x=S.OHm) 




EsUnu(|ed valu^ point A (x=3.0Kni) 




Figure 5; True and estimated reflection coef- 
ficient plotted against reflection angle at 
point A and B. 



Having generated the input data to the experiment, all the 
ray data needed for the amplitude correction are com- 
puted by ray tracing by wavefront construction which is 
fast and reliable also for 3D cases. 101 stacking sticks 
with length Lp = 250m and a spacing of 100 m are 
placed along me target reflector. The stacking of syn- 
thetic pulses from each of the 10 .offset intervals as 
described by Eq. (2) gives one i^synthC^) ^^^^ 
interval. Following the method description above, w.e end 
up with ten estimates of the reflection coefficient 91^ (<P') 
for each of the 101 stacking sticks. 
Each of the 10 estimates c.orresponds to a specific reflec- 
tion angle. In Figures, Sl*^(q)*) in two x-positions (A and 
B indicated in Figure 6) and 10 reflection angles are plot- 
ted. The plot shows excellent correspondence between 
true and estimated values. 

Comparisons of true and estimated reflection coefficients 
may also be done for constant reflection angle and vari- 
able X as in Figure 6 where for the near offset 
(25-300m) is plotted along with the true, normal inci- 
dence reflection coefficient. Also here, the method gives 
good approximations, but, as expected, the resolution 
varies a great deal. Note how well the two weight func- 
tions at the top of Figure 6 (W^ and at point C and D 
respectively) describe the "smearing" of the true reflec- 
tion coefficients. 
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Figure 6.* True and estimated reflection coef- 
ficient of normal incidence rays. Two of the 
weight functions giving the resolution in each 
point are shown at the top fWi arul Wa^. 
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